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ABSTRACT 


The heavy dopiiig phenomenons of handgap narrowing 
and Auger recombination have been numerically modeled. 

Their individual and relative importance in detemining 
the current gain of an n-p-n power transistor have been 
studied* It is shown that in the devices with emitter 
junction depths of 2{am or less bandgap narrowing dominates 
all other mechanisms . SRH recombination has the most 
effect on current gain in the devices with emitter junction 
depths of or aiore • 



CHAPTER 1 


IHTEODUCTIOE 

Tile opera.tion of a junction transistor depends on 
the introduction and transport of minority carriers across 
the base region of the transistor. In the design of an 
effipient transistor, it is necessary to contrive that the 
emitter current is predominantly a current of minority 
carriers emitted into the base rather than a current of 
majority carriers withdrawn from the base, Por this require- 
ment to be met the emitter side of the emitter ba.se junction 
must be more heavily doped than the base side. 

In most silico_n devices, the emitter doping levels 
in excess of tO are found. At such high doping 

levels, the basic parameters such as the density of states, 
the width of the forbidden energy gap, the minoritjj' carrier 
mobility and the minority carrier lifetime strongly differ 
from their values in a lowly doped crystal. The bandgap of 
silicon is reduced, the electron and hole mobilities decrease 
and the aiinority carrier lifetime becomes lower at heavy do- 
ping levels* 

These heavy doping effects have to be taken into 
account in the modeling of transistor* The electron and 
hole mobilities in silicon as a function of dopant concen- 
tration and temperature are important parameters for device 
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design and analysis, The importance of the minority carri- 
er lifetime is obvious , It influences dc as well as 
transient perfoimance* The impact of bandgap narrowing 
effects on device performance is nearly as obvious. The 
bandgap shrinkage and the deformation of band sti-ucture 
results in an increase in the pn product with increasing 
impurity concent ration, Zeopiiig 3.n mind that most electrical 
quantities of interest, such as currents and minorit;r ca- 
rrier storage, are proportional to this product one under- 
stands the impact of these effects on the device performance. 
Heavy emitter doping has a detrinxental effect on many tran- 
sistor parameters. In addition to decreasing the transistor 
gain and increasing the temperature dependence of gain, 
it also results in a higher noise figure* The excess noise 
generally attributed to surface effects is substantially 
larger for heavily doped emitters, 

1 ,1 Literature Survey : 

Mertenset have considered the heavy doping 

effects by modifying the device transport equations. They 
have calculated the transistor current gain after solving 
numerically these modified transport equations, Slotboom 
et al have given an emperical fomula for bandgap shrin- 
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Ijage, which, can he used fo: concentrations in the range of 
4xl0^^ - 2*5 Lanyon et have derived 

an analytical expression for handgap narrowing which may he 
used for modeling of the handgap shrinkage* Ihe expression 
is in good agreement with experimental results in the doping 
range of 5 x to t ,5 x cm”^* Possin et al have 

suggested some modification in the classical device equations 
for including the handgap narrowing* They suggest using an 
effective electric field term in the equations. G-aur et al^^^ 
have modeled the heavy doping effects hy including Auger 
recomhination and correction 'co intrinsic carrier concentra- 
tion* The results of measurements of handgap narrowing done 
hy a numhor of investigators is shown in Pig, 1* 

[71 

Shihih et al*- have shown that Auger recomhination 
alone cannot explain the common emitter current gain in 
bipolar transistors* Other physical mechanismsin addition 
to Auger recombination have to he considered to account for 
the observed values of current gain. The Auger coefficients 
calculated hy various authors are tabulated in Table t * 
ITumeriftal treatment of Auger recoEihination has been done 
by ~ Eecently Polsky have numerically 

simulated the bipolar semiconductor devices taking into 
account heavy doping effects and Peimi statistics. 
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t .2 Preview : 

The objective of the present work is to siinulate the 
heavy doping effects of bandgap narrowing and Auger recombi- 
nation in order to deteimine their individual and relative 
influence on the current ga-in of ann-p-n power transistor* 
fhe basic device equations in one djjnension with required 
modification for inclusion of bandgap narrowing and Auger 
recombination effects, have been solved for the transistor 
and its current gain for the given impurity profile has been 
determined in the presence of different combinations of 
these effects. Chapter two presents the physical basis 
for the modeling. Chapter three discusses first the general 
philosophy of modeling and va.rious approaches to the modeling 
of bipolar devices. Iho approach followed for the aiodeling 
of n-p-n power transistor is tiien described. Chapter four 
presents the rosults obtainod and the discussion of the 
results. Plots of current gain versus emitter junction 
depths and current gain versus collector currents are presen- 
ted. Chapter five summarizes the main rosults of the work. 
Ihe computer programmes written in Fortran have been included 
in the Appendix. 



CHAPTER 2 


THSORATICAX BACEG-ROUHE 


2*1 Energygap : 

At heavy doping levels, the numher of naohile carriers 
is very large and they are q.uantim\ mechanically indistinguish- 
able i.e. one particle cannot he distinguished from the 
other* The energy distribution of carriers in such condi- 
tions obey the Eeimi-Dirac statistics rather than Maxwoll- 
Boltzman statistics which is obeyed at lower doping levels 
when the particles are distixiguishable * 

Since there is a large number of carriers in the 
crystal, the separation between them is very small and the 
interaction between carriers, and carriers and ionized ions 
becomes prominent. The majority electrons (for n-type Si) 
screen the donor ions through Columb interaction and reduce 
the ionization energy of the donor* This results in donor 
level moving up toward and ultimately into the conduction 
band* Similarly, when majority electrons interact with mino- 
rity holes, the hole potential energy is reduced and the 
valence band edge moves toward conduction band* The electron- 
electron intera.ction is due to their spin motion and there is 
an exchange energy associated with this motion. The effect 
of this interaction on energy band is that the conduction band 
ddge moves downward. 
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Ilie distribution of impurity atoms in the crystal at 
heavy dopings is not unifoim, The crystal is disordered and 
the distribution of donor ions is random* This random dis~ 
tribution causes fluctuations in local electrostatic potential 
and variations in the energy of the electronic states in both 
the conduction and valence band results. These variations 
in energy cause a spatially dependent distortion of quantum 
density of states, Concommittantly the statistical average 
over the entire lattice, of the density of states, -which 
defines the macroscopic properties of the semiconductor, 
sho-WBtailing into the energy gap, of both the conduction and 
valence band density of states* The bands are not parabolic 
near the extrema in such cases and the energy gap is effect!-, 
vely reduced* 

The discrete donor levels at low dopings levels, 
spread to form a band at high doping levels , This is because 
the spacing between the impurity atoms is small. This allows 
the electron wave function at the bound states to overlap, 

!Ehe random distribution of the donors also results in vari- 
ation in the donor energy level and formation of impurity 
band tails. The effect of the impurity band formation on 
band structure is negligible since the impurity bands are 
narrow, there is only a small aumber of electrons in them, 
and the mobility of these electrons is small. 
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Th) msasuromants of 'bandgap narrowing have yielded 
two typos of tandgaps. The mcasureDent of optical ahsorp- 
tioa give tho ’optical’ energy?" gap. Tlio change in this bandgap 
is due to many body effects such as tho interaction of elec- 
trons with electrons- and with holes. Tho meesureiBsnts of 
minority carrier concentrations 3rield the ’electrical’ energy 
gap. The change in thif^ go.p is due to both many body effects 
and band tailing effects. This change in ’electrical’ energy 
gap do fin -s the effective intrinsic carrier density n^^. 

2.2 Auger Pu.^ combina tion : 

The band-to-banc Auger r combination is a throe carrier 
effect, in this px^'ocoss, tho recombining o-h pair gives its 
energy either to an electroiu in the coxiduction band, (called 
e~e~h process) or to a hole in tho valence band (called e-h-h 
process) . The Auger recombination predicts lifotirio inversely 
proportional to the square of tho majority carrier concen- 
tration in degenerate as well as in non-degenerate semicon- 
ductors . 



CHAPTER 3 


DEVICE I-TODEIINC- 

Modeling is tlie art of chars-cterizo-tion of the beha- 
viour of a physical process, Ihs device ffodol attempts to 
describe the internal and the terminal electrical behaviour 
of the device. There are two broad categories into v/hich 
modeling can be devided* Both these approaches are comnionly 
used in device modeling. The parameteis of an empirical 
model are generally obtained through a. study of the e^cperi- 
mentslly m *asured behaviour of the device as Yiewed from its 
temlnals . This involves application of curve fitting tech- 
niques to obtain functional relationships between the teraiihal 
quantities of interest. The physical model, on the other 
hand, is based on an analysis of the basic internal physical 
mecliaxiism of the device* It enables an ixisight into the 
working of the device and gives an excellent qualitative 
picture of the effects of changes ixi material and structural 
pa.rameters on its teimina.l, electrical heha,viour. 

Models have also been developed which make use of a 
combination of the two approaches aientioned above. The 
physics of the device is used as a guide to its terminal 
description. As far as the device engineer is concerned, 
it is the physically oriented model that is more important. 
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This is heca-use only this type of model enables the opti- 
mization-designing the device so as to improve some specific 
aspect of circuit performance to be performed. 

3.1 i s t o r_ I'T 0 d elin g ; 

The following throe basic approaches ai’e in use for 
modeling a tra,nsistor,. 

(a) The cla.ssical ana,lytic approach : The models developed 

using this approach are simple, can be ans.lytically maaupula- 

ted, and are directly z'elatod to the physics of the device. 

However, they do not accurately reprosont the behaviour of 

the device under all conditions of operation. An e^cample of 

this type of approach is the w8l3--lmown Schockley iDodel of 

( 

a transistor, 

(b) The total compuh-.r a pr-ro ach s This ap'oroach is based on 
a rigorous soluti.on of th"; basic semiconductor equations. 

The results obtained are highly a-ccur.ate, 

(c) The ^regimi^ a;^roxiiLmt_^on method : This approach is 
based on the division of the transistor structure into well- 
dofinod physical regions. Different a,pproximation are made 
on the basic semiconductor equations in the different regions 
and these equations are then solved either anal 3 rtically or 
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on the computer. The self-consistency of the approximation 
is checked from point to point. 

The present n’odeling follows the tota], computter 
approach. Th' hasic semiconductor equa.tionsare solved hy 
finite difference method after including dto edifications due 
to heav 3 ^ eloping concentrations . The handgap nar:x5wing is 
modeled as correction to the intrinsic carrier concentration 
and the modeling of Auger recombination is done by adding 
an expression of Auger recombination i-ate to the expression 
of SES recombina.tion« Por mobilities expressions are used 
for its dependence upon carrier concentration ^■,nd electric 
field. 


3*2 Prob lem Definition ; 

The problem is to obtain current gain of an n-p-n 
power transistor with varying emitter junction depths by 
solving a. set of three basic semiconductor device equations, 
introduced in the next section. The following information 
about the device is given* geometry, doping profile, mobili- 
ties as function of electric fiild and impurity density, 
genero.tion recombination laws including Auger recombination, 
intrinsic carrier concentratioja corrected for bandgap na- 
rrowing and the appropriate boundary conditions . 
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Sa-sic Equations : 


Tlio sepiicoxiducoor tra-iisj^or-'t "jq.up/bioia'a which, have to 
he solved for thj steady state behaviour under isothemial 


conditions a 

.re 
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Equat 

ion 

(1) is 

tho Poisson’s equo-tioa, 

relating tho 


divergence of the electric flux to thj oluctrostatic charge 
due to nohile holes (p) and electrons (n) and tho ionized 


impurity atom density (N = - Np . Tho 

(3) are the continuity equations for the 
under st-^ady state conditions. Equations 


E q uat i 0 as ( 2 ) and 
hol.iS ?nid electrons 
(4) and (5) are the 


current density relations composed of drift (duo to tho 
eloctHc field) and diffusion (duo to car'tLer concentration 
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gradient) tiims, Boltzaian approximation of Pemii sto.tifstics 
■between the potentin.1 and carrici* donsitiwS have 'been assumed* 


{ ^ f ) 

n = n^ exp (a 


(7) 


P 


n. 


exp (q 




( 8 ) 


Quantities co^ and 9^ in 'Equations (7) and (8) represent hole 
o-nd electron quasi-Foimi potentials. The Einstein relation 
betw >en tlio carrier mobilities anid the diffusion coefficients 
are assumed to bo valid, although this cannot be clone for 
m.ajority car-’i 'rs in degenernte semiconductors. The depen- 
dence of caricicr mobi3.it? es on electric fielC ana impurity 
atom consity is mod.eled a,ccording to the following formula. 
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Tile SlioclCLey— Road— Hall recomliiiation law is given as follows 


Rl 


P-n “ 

■■'no'^iTT + tpirsT^' 


(11) 


Tile equations for tlie Augar r-ocomTaination aniC higli field 
generotioh rate are the fo3-lowing 


R2 = (A^n + ApP)(pn - np 


G 


n + P 4^ 


P P" 


(12) 


(t5) 


Tlio va"’ ues of v;.riouB constants entering tlie eqviatlons (12) 

LH'4 

and (I 3 ) arc given in Triolo 2, 

Tlio tra-nspoi’t oquatioxis (1) ~ (6) arc normalized into 
dlmenslonlee® forms according to noima,lization factors of 
Table 3* Tlie iioimalized set of equation is reduced to three 
coupled equations in variables 'if n and p* 




n-p-ll 


(U) 
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4 Yp‘)) = R - G (15) 

-7a)) = R - C- (16) 

Bqua'cions (7) and (8) in their noimaliZsi'd form reduce to 

a = exp (^-cp^) ( 17 ) 

p = exp (cpp-4') (l8) 

Here two vazlables and are explj citly introduced 

$p = exp (cpp) ( 19 ) 

£5^ = exp (--cp^) (20) 

These expressions are substituted in the equations (l4), (15) 
and (16) to obtain a n-?vr set of eouatioiis in the -variables 7, 
ajp and 

7^f = exp(1?) - 5^ oxp(~W) ~ R (21) 

ir,(pp exp(~4?) T ip) = R - G (22) 

7.(p^ exp(4f) 7 i^) =R-G (23) 


These three coupled equations have been solved by the method 
outlined nvjxt. 
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^ cretj.zati on ; 

Tlio Poisson’s ecj.uation (21) is nonlinear in Y. This 
is linearized "by assuining 

6 (24) 

wliero 4 is the potential value of the previous iteration and 
6 is th:. correction toim* ¥e have 

+7^6 = oxj){'$ + 6) ~ exp {-'^ 6 ) - ¥ 

= i exp( f)[ 1 + 64-0(6^)]-$ exp(-.W)[U6+0( 6^)]~N 

Jr 

= exp(^) -ipOxp( ~f) -N4 6( $^oxp( ’5’)45p3xp( ~7) ) 

. . Y^5-6(S„exp(7)+i oxp(~w)=_>y 2^+S oxp(W)-S exp(-^)~iT 

iX T) J.X ^ 

(25) 

Thus neglecting the second and higli3r orc'or toia^s 
equation (25) is considered .ns a linear differoniial equation 
for 6, 

The discretization process is achieved with the 
selection of a unifonn step distrihution. The finite diff- 
erence scheme used for reducing the Poisson's equation to 

the equivalent lin.ar algebraic equation replaces the first 

"bh, 

derivative at the i step by the formula 



t6 


du 

dx 





2h 


u ^ 
..JL-J, 


and ’the second cTorivati’^e "by 


d^u 



i 


Vi \ 

2h^ 



( 26) 


(27) 


where <^^41 » "values of the function 

ij *H *ij *h. "ij 

is at poin-isi+1 , i~l , i '' distsaice between two 

successive points, respectively. Using these formulas, the 
difference equation at the point ’i» [ilg. 2] is writ'uon as 

[2 + h^ (i^^(i) exp(^f(i)' i* 35^(i) oxp(~’^(i))] 6^ 

= + exp 

•“Cpp(i)exp(~ f (i) ) - U(i)) (28) 

The continuity equations are replaced by the diff3X'enc3_ equa- 
tion wh 3 .ch are obtained by the box integration method. 

The difference equation corresponding to the equation (22) 
is given as follows 

( a( i+i) +a( i-i) ) l!Sp( i ) =a( ii-i) i^( i-h 1 ) 4 a( i-i) i-1 ) 

whe I’e a = Pp exp ( - 'F ) 

The value of *a’ at the point midway between two 
successive step points is approximated from the taiown 
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■values of 'p.,’ af tluso fwo T^oia'fcs . Tlio differoace ec[uatioj:i 
for the rola.tion (23) is similar to (29), hut with a=ji^ozp(^). 


3 *5 Sol utjo a Procodurj : 

Ihe differonce oquations derived in t-he last &?ction 
are applied to evoi'y point of the step distribution -sj-hich 
give ri '9 to a S3''stem of simultaneous algebraic equation. 
These equations can be written in a. mat rise forni as follows 


AJ = Q 

This matrin equation represents a syste^'' of 
where each component of the unkiiown vector 
only to adj icent components. 


( 30 ) 

lin-'ar equations 
’X’ is related 


The matriz thus foimiod is of triple -diagonal form. The matrix 
equation is solved directly by a, method known as Grout’s 
method- It is a numerical procedure for forward elimination 
followed by backward subs'fcitution. The recursion fom‘'ulas 
for tho technique are as follows : 





= I’m + \ ^-1 

m=n-l , n-2 , , . 

... t 


( 32 ) 
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whe ro 


la = u 

m ^ b la . » 

m m m-l 


Q - b i 

_ in m 

” “b +■' \ V-1 ’■ 



(53) 


Tile seciuence of solution for ^ , cp and cd^ is similar 

P 

to the one descrj.bed in reference [1^] and is cEiled G-umniel’s 
algoritliD\ Pig. (3). 


1 , An initial guess to T is made fx-om whicii initial 
valuosof cpp and cp^ are calculatad for tlio given anjlied vol- 
ta^3e. 

2« Poisson’s equation is solved and new value of T 
is calculo.ted# 

3m Continuity equa^tions are solved to obtain now values 
of 9 p and 

4« A test of accuracy is done. If required accuracy 
is reached the iteration is stopped othervilse go back to 
step 2, 

The cui'rent gains aro calculo.ted as the ratios of 
collector current and bu.se teminol cui'i’eiit , The effect of 
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various pliysical B\ocli£.nisms with. vo.rying 
incorporutf?d by doing the ca3-CUlation on 
lin ''S w5.th cinitters of diffGr''nt sizes, 
computations the total emitter charge is 
emitter depth is varied. 


emitt:-::? size is 
the above mentioned 
For this type of 
kept constant as the 



CHAPTER 4 


iGSELTS AHD DISCUSSION 


4*^ Res ults 


The calculations wore pr3rfo'miGd for a poorer transis — 
tor structure which has the following douhlo Gaussian doping 
prof 5.1 e • 


l\r(x) = e 




Oe' 




where the distance x is meosured froDi the surface 
transistor. The structural and profile po,.rariieters 
transistors studied are given in Table [4]#. Rig* 
plot of the doping profile* 


(34) 

of th^- 
of the 
[4] is a . 


The mate.rial paranetors of the transistor are given 
in Table [ 5 ] . 

The dcteimination of the emitter surface concentra- 
tion r-squired for oqual emitter chcirg.-i for vai-yi^ig emitter 
junction depths was done using pTOgrame (2), 


The results of the c aloud. ations performed for the 
detoOTination of relotive importance of various physical 
phenoDienons by varying the emitter depth are tabulated in 
Table [6] amd sijjcma om Rig. (5)* 


The resiilts of the computations performed for the 
doterriiination of curi*3nt ga5_n variation with colic cito&Al^l^^nt 

'I LMMiAm 

^ •’S- A, 
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in oho pr^^Reiice oic dj plij^'Gical niechanisiiiB are reported 

in Table 7 Ji,nd plotted in Pj.g. (6), 

4*2 DiscuGBion 

Fig .(5) shows the depondenco o;C current gain on 
enittor junction depths in the presence of different combi- 
nations of physicol processes. The current level in all 

2 

cases is 8»2A/cm , First, the case of 4pJ® device is consi- 
dered since it is typical of the overall results. The current 
gain in th'; situ'viion when a3_l the m "clianisins are contribu- 
tiijg, is just bcloi'r lon, if b'^^ndgap narrowing is removed 
now the cur'’ont gain tncreo/.--'S to I 40 , showing no significant 
infd.uenco of bandgap narrowing for such a, device. If bandgap 
nar 'owing is rostorod and emitter SEH traps are removed, the 
current goa* n increase to 425* The removal of Auger x-ecombi- 
ziation in com.biu-’tion of an^^ of the other physical phenomenon, 
produces only a litti.e furthr.r change in the current gain. 

Thus this shows that SEII recombination is tho dominant mecha- 
nism influjneing tho current gain of the devices with emitter 
juxiction depth in the range of 4 ihn* As the junction depth 
is further increased, tho degree of this effect is also 
increased. 

In the ca,se of the devices with emitter depths of 
2p.jri or loss, tho direction of these dependencies is s'ame 
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but tho rol.rbive magnituur.s change. Kio banag:ap nax-:’Owing 
produces larger chraig -js in tb© cu^-ront gain on removal, The 
rJDIi recombtaabion has littlo effect on currxent in this range. 

The dop.sndence of cui'rent gain on emitter ^junction 
doprh in the prosonce of all m-cliajiisire can now bo explained* 
As the emitter depth is incroasxd from 1 to 2 pm the 
bandgap n.ar..’Owing effect weakens s,nd therefore curronb .gain 
incroo-Ses • In the devices of en'^itter depths gi-^ater than 
2 pm, current gain decreases with the depth since the effec- 
tive emitter charge is lovrer becausv) the net impurity density 
is giMdually decreasing in the active emitter region near 
the bese region. 


Thgures (6) shows th.-, vnxd.ation of current ga-in with 
increasing coll-;ctor current I.VvjI. Curve (a) shows the same 
for th<j case when 'll! mecbanifmis are oper-'tingj Ourye(b), 
when Auger r cor'' bination has been removed Ourvo (c), when 
bandgap narrowing has b ) ni removed and Curvo ( d) , whom 
bandgap narrowing and Augor re combj nation both have not been 
concidorod in the analysis. It is clear from the fig;> that 
Augor r',com,bi' nation does not have much influonc.^ on tie 
current gain. On removing bandgap narrowing 'fz-om the analysis, 
the resulting current gain increases, Removing Auger recom- 
bination alone or along^rith bandgap narrowing produces little 
further change in the current gain. 
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Tho impo 2 '’''cai;ico of SlflT ro combination in date lining 
curi'ent gain can bo argued as folloi’^ » If tlaere is no ro~ 
combi mat iou, baadgap narrowing has virtually no effect on 
cur'-’cnL gain with inci-’easing inj action level* Bandgap 

nar 'owing enhances tho rocombination rate if recombination 
is occuring and this reduces th;; current ga.in. As is obvious 
from th'j oarlisii- d.iscus>sion, this is tho main effect cf band- 
gaT' no.rrOvrinG in devices with emitter junction depths grea'cer 
them 4 p.m * Augor recombination increases with injection 
lev j1 thus reducing bhe current gain or cit best leaving iu 
unchoniged* iligb injection level effects also tend to reduce 
the current gain. Contrary to all this, SRH recombination 
deer. xrith inj jction I3V0I and thus is tho oni.y mechanism 

th-.t can e:qaain th: observed initial rise in the current 
gain with injection livol. Bandgap narrowing can affect the _ 
i.>combinat-iou rate but alone cannot explain the observations. 
Those argumeutf? are confix^ied by the necessity of adjusting 
the emitter lifotime at different injection level to correctly 

predict the cur''’ont gain. 
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TABLE 3 


V. -ilaLlo 


Normalization Table 

Nomalizod Factor of normn-lization 

(quantity 


Llcctro.'itati c Pot nitisl 

'^ua.o i -!■'(, i n i I o v ol o 
Volt. 

p OG i 'b :i 0 n coo rcli n \t ■; d 

Ch'-.r, :f. Cone ;ntr.-'.tionrj 

Ca;'"i r L? ' "uiiion 
Co* ffici.jntrs 

0. ir*:.’!'! r m nlil It i jn 

Cur :n nt doiirs -i.t n ; iO 

C!’c n. ) r • ' t i 0 j 1 ro c ombina- 
tioii rate 

Carrier lifetimes 


V 

U'i. 

^n^t ' ^p'^t 
Pn » l^p 

R, Rl,E2,G 

'"no» 


kT^ 

V “r = *026 ^ 


V. 


V^ 
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S=jalx7 = 


1 Pi 

n^=1 .5x1 0 ^ cm 
L^ = 1 cm^/sec 

2 

--• = 39.06 cm /sec V 

H 

n p 
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D c nf 
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Fig.1 Bandgap narrowing measurements 
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comusioij 


It is sho^rji that tho hojidgap iiarroifiag £,ncl SEH 
rocomhination effects are tho dori''inaiit nochaiiisn^ in influen- 
cing th'- curvemt gain in transistors, Bandgr.p narron‘5.ng 
effects dominate nl?. others in tho devices with sh'^llow 
emitters of junction depths of 2 pm or less, The SPJi recom- 
hj.na.tion has the most influence on t’le current gaj.n of the 
devices with ero' bters of 4 pun or greater. Auger recomhina- 
tion does not appear to ho an important mecha-nism at1?he low 
and modnr.ate injection le^ols. Ihis is dominated hy either 
han4SS'P narro'ving or 3HH recomhiaation and ca^nnot explain 
tho ohseir/od currnit g.'.ixi in t]ie ahsence of the above two 
mechanisms, Ihcne conclusionc havo b-exi arrived at by 
choosiiig specific values for th-' pararnoters governi v, the 
physical mecha.dsms,for example, b<':uad£;ap nar-oowing, Iho 
specific choice, however will only effect tho results for 
very small emitter depths oir»ce x’ecombination prevents the 
minority carrier from reacliiiig th^ portion of vhe ei'^i'cter 
doped at t^iis level for the deeper jmjtter (2‘-4p-m). 
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Appen dix 

Solves tile device equations 
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Pk 3 G ? i^v f , io ,3 

FhJS PROGRAM DFiTF.RMINES THE COMMON £:v!r:rrER CURRE^jI GAI.N FOR 
POftER TRAMSlbTOR AFTER SOLiVLNG TRS'. 3ASi:. SEMlCOMDUCl'DR 
.RrtMRPJRi* E.pATID^S USIiMG^ FINITE DIFFERENCE: METHOD, GUMMED ' S 
Rt.GJkJIHM Af^D CRGiJT'S FACTORISATION MS, rH3D. 

VAKIABDE DECLARATION 



, E(9D0), P( 900) ,B( 900) ,C(9D0) , DODO) ,3(900} ,DELC9D0),DSLMA,DE 
iLHN,EFC9DO),FU9o6),R2(95o),v!UPC900),MJ.»i(90a5,G(900) , JN(90!)/*JP 
4(900 ), hfe,e:psj,x(9o6) 

INTEGER' i,J,iri,Z 
DUTPUr 


DPeSi(DMir=30,DEVICE='DSK',FILE='PP.T ') 

NORMALISATION FACTORS 

3 / 1 = 0,026 

L‘D=34,e-4 

.41 = 1. 5E1& 

D 0 = I .', 

Mils 38. -65 
J0=7,E-7 
10=34, €-4 
AA = LO * 4'5 
APPLIED' BIAS 
VA=0.0 

ITERATION STOPPING CRITERIA 

inM = 200 

EPSi=0,lE-O4 

CONSIDERATION OF DIFFERENT MECHANIS'MS' 
Z=.l FOR NOMINAL CASE! 

Z=2 FOR AR REMOVED 
Z=3 FOR BGN REMOVED 
Z=4 FOR.BGN+AR REMOVED 
2=5 FOR EMITTER SRH REMOVED 
, 2=1 

DOPING PROFILE PARAMETERS 

NE=(7,E+19)/NI 

MB=(2.E4-17)/NI 

N:,= (2.€.H)/Ni 

&ETArI,-2 

ALPHA=ALOG(NE/(NB>t=EKPC-BSTA)-.NC)) 

KE=C4.75E-4)/LD 

JC&=a2.25E-4)/LD 

XC3=(90.E-4)/LD 

STEP DlSTRiaUTION 



P ft ji£i 


'■ii = (u.u>4)/La 

ri - 4 1 I' H ^ 

^ - X Z/ H i + i 
DD lu J=l,f. 

XCl)=ItHJ 

.lf^*^^i^^4ISA'J‘XD^i OF 'JhRlABUES 
JJ 7'J 

oei,ci)=o,5 

PHI?(I)=i.O 

pHi\cn=i,o 

C'JNriNU’E 

CSiLL^ D3PING(M,NE,NB,fjCf 

«IEFr(N,VT,L^i,DOP,E^’] 

"Ai.r. 'STTKiTrw c'iiT hna at \ 


W »wtJ' JRi Vf Cj 

Cs^LL. D3PIMG(M,N'E,NBtfJCilLPHft.BErft,X£> )t,[>3P;) 
CftGL NIEFFtN, VT.L^i, d5p,EVI,) 
CftLL^SIIMT(N,ENf,6op;SI) 

call' CCEFFl (^,H2,VA,PHIP,PH[I.N,SI:,DD?^, Bi, :;,D> Q) 

ZAiy crout(n,6.c;o,5.del) 

CALL' DL«^(N,SI,DBLj,6eLMX,DEL'4N) 
IFOELMX.LT.0.1E-04JG0 rC 420 
IF( III .EQ.ITIN) GO rO 430 
03 410 J=!,N 
SIC JO*SICJ)+OEL(J) 

CDNTIMUE 

iTi=.iri+i 

G3 rO 405 

WftirE(30,9000)ITl;G3 TO 440 

WRirE(30,910)Il’lM 

03 450 I=l/N 

sici)=sici5+del(I) 

coNiiNyi; 

CALL* PHf(N,PHIp.PHIfl,III,SI,P,E:,PHP,PH>l) 

CALL* ELFD(f»,Hl,VA,SI, EF) 

:ALLi RSPH(f«,Hl,XE,Z,P,ErMI,Rl) 


:ALLt ELFD(f»,Hl, , 

ZALLi RSPH(f«,Hl,XE,Z,P,ErMI,Rl) 

:alli ragr(n,p,e.ni,r2) 

SaLL* MCBLTYCN.MO.DOP, IFj-HUP.^UM) 
:ALLi flFG(N/P/E;M(jp.^UN,ErfG) 

:ALLi C3EFr2(N,H2,vLz,MUPeSI.,Rl, 
:alli C'Routcm,b,c,d.3,phi.p) 


G) 

SI.*Rl*R2,G.*3i,Ci,D,3) 


200 STOP. 

9000 FOR'^ArttOX, 'POT 

919 FORVIATCIOX, 'NO 

END 

;i , THIS' SUBROUTINE 


•a & r r 4 ; « £ n 4. f t ri # 4 * # « w r ■£ 1 , .# n 4 . # r\ 4 # a»# h j 

3EFFl(S,(i2;SL£,!15s,SI,Rl ,R2, 3,3), 

? “W A n . ' ■ 


u r f. ij P f 1 r n 4 f ft f n r /■ I If* r\ X A .? V f / XI r & 

POTENTIAL SfALUES CDMVERGEO- AFTER ' , 1 3 , ' I TERATI 3NS' ) 
NO REOUIRED CDN^^ERGE^ICE. AFrER' , 13 ITERATION ' ) ’ 

INE DETERMINES OOP! NG PROFILE: 



^ 1 / 
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3 D 

20 

10 
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10 


a'lt^RJi/riNE DUPIflG(N,NE.,M3,M:;.ft[iPHAt 3,e,rA, Kg:, X,DDP) 

01.'a£:iS13w D0P(N),XC^), £1 (9003 , E2(90D3 
HKaL 

0 3 1, w I = 1 , 0 

£1 ( U = &.LPHA>^C (XCD/XEO^^a) 

IF £1(13. GE. 20) GO TO 20 
&2( I.3=8F.rA’»‘( CX(I)/XE)’i'^2) 
OOP(£)=^JE*EXP(-EUI))ftC-N8»=EKP(-E2Cl)3 
w J f 3 1 0 

£2n.J = 3ei’AnCX(I)/KE)»= + 2) 

iF( E2( U .GE.20)G0 TO 30- 

i)3P(I) =.-fiB=l'EXP(»E2( I) )+<M:i 

G3 f3 10 

03P(I3=IC 

ciNTi^ioe 

RETOR?^ 

£’?iD ' 

rms SUBROUTINE DETER^INSS 3.GN OORRSOTION FACTOR IN Ni: 
.SUBROUTINE Nl EFF ( N, IfT, Z , SI , OOP , ENI ) 

DIMENSION DOPCN) ,ENlCb,FC§00),AD3?C9003 

REAL.- NI 

03 10 1 = 1 ,N 

IFCZ.EO.SIGO TO 30 

IF(.2,e5.4)GO to 30 

XL=C1.E17)/NI 

ADOPCr)=ABS(OOPCI)) 

IF( AOOPCD.GE.XU) GO TO 20 
EMI(I)=1. 

G3 TO 10 

FCI) = AUO3(A00pa)/XLD 

£NI ( 13 sEXP CC9.E-3)*CFCI)lSQRrCF( !)<'£( 13 4-0.5) 3 /VT) 

CDNTINUE 

RETURN 

END 

THIS SUBROUTINE DETERMINSB INITIAb /ALOE. OF ELECTROSTATIC 

potential CSI) 

SUBROUTINE SIINT ( N , ENI . OOP, SI ) 

DIMENSION ENKM) , 06 pCS 3 ,SICN) 

03 10 1 = 1, N 

IF(DOP(I).GT.O,03GO TO 20 
SrClJs-ALOGC-OOPd)) 

GO ro 10 

SIC I.) = AL0G(00PCI)> 

CONTINUE 

RETURN 

END ■ ■ 

THIS SUBROUTINE DETERMINES CDEFFIC I EN T^ MATRIX FOR VARIABLE! 
SUBROUTINE COEFFl C N , H2 , V A, PHIP , PHI N , S L, OOP, 8 , C , 0, Q) 
OIMSHSION PHIPCM) ,PHIM(M3 ,Si:(N3 ,OOPCN),B(N) ,CCN3 ,DCN3 ,Q£N) 


DttL 



§ 0 


23 
33 
I 3 


ID 


43 

23 

33 

13 

:i 


3 D 1 0 I = 3 , N ■ 
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33 ro ID 
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THIS SUBROUTIME DETERHIMSB SRH. RECOMa£..VAriO!4 RATE 
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REAL* fll 





j 

Si 

1 D 


23 


13 


pi iu l=i,Jl 

i F ( 4. 


■'(4.f:J.5 3(;D TO &0 
!u=.5,8E-6 


f Pi# "5 » 8 

iji D F J 5 3 
i'NU='i .5fc:”.4 
i'Pu=.i .5E-<5 

gu f.J=CPCI)*ECl)-f'ir#NI)/C r?i3*(p(l)f^l) trpD^CECD+^I)) 

», .) I’ i. U E 
Jll=JUl 
03 2'0 i = jn,N 
r '!U=I,'5E-'^ 

r?U=.1.5E-4 

kU I.) = CPCU’l‘EC l)-Nl^NIO/CTNO*(PCI)f ^I)frPD*(ECI) + NI>) 

CJ^riNJE < 

RET’JRM 

ENfD 

THIS SUBRQUTIWE DETERs4II^SB 4.U3ER RSCD'liBISAIIDM RATE 
SUBROUriNE RAGRCN,P,E.NI,R2) 

DIMENSION P(N),ECN),R2(N) 

REAL MI. 

AM=3-.1E.-19 

AP=2'.2E-19 , ■ 

03 10 I.= 1,M 

R2CU = ( AM»E(I)tAP^P(I3)’f=CPCr,)<‘ECI) 

C3.MriMU.E 

RETORM 

EMD 

THIS SOflROUTINE DETERMIMSS '^OBILiITI E.S 
SUBROOriNE MOBLTYCSrMUrDOPfSFfMUP, MJ.M) 

DIMEi^SION DOPCN) ,EF(N) .P'4UiC9DO) ,E'4j;i(900) rFlC900) , 
lF2C900),F3(900).Gl(9o6),G2(900)rG3-(93 0) 

REALi '<LI,MUP(N),MONCM) 

D3 10 1=1. N 

FUI,) = COOP(I)/(CDOP(I3/8l.)f2.7E6)) 

F2a) = ((EFCl5/7,9E2)**2)/C(EF(I)/7.9E2!)4-,2l) 

F3CI)=C CEF(I)/3,2E33**2) 

PMO! i I) si . +F1 ( T 3 + F2 f 1 3 +F3 { I ) 

P«^Ulh 4: + (D0PCI)/ C&OPCh/81)t2,7SS))H(CSFCI)/7.9E2)*»2‘)/.UEF 
HI )/7.9Eb+0.2ipi-((EF(I)/3.2E3)* <<20 
WPCI)=.C(PMUHl))*=*'No,Sn<*'f80/39 ^ 

ceF{li/4l5E2)Ji^)/C?EFCif /i.«S20 4'l .'l93 


G3CI0=CCEF(lj/9.6E2)**2)^ 
EMUia)=l.+Gl(l5 + G2(l)+G3(I) 
EMUt a)=l. + COOP(I)/(CDOPCl)/ 


_. + COOP(I)/(CDOPCl)/350)+2.'S:5)) + ({CEF(I)/:4.5E2)»»2)/.HEF 

HI)/4.5E2}+H15))t(C|Fa:)/9j6E2)*«20 
MUMC I) S.CCEMU1 (!)}♦* (-0.5) )^l400/39 



^ A ^ . 


3 ») 


i J 
1 J 


4 3 


50 

20 

50 

10 


cowfi.siue 

t(!.%rj.HNi ' 

il.'iO ■ ■ 

i‘(iJ3 SJBRDUTINE DErSRmMSB -l.1 3H FISLO' 3S:!^£RATION RillE 

cjilbROlU'fNE hFGC^^,P,E.MUP,MUN,EF',G) 

0/f4E,NS1.3i^ PCi.) ,E(S) ,mOP(M) , ^U.MCM),S?C >i) f 3-(M) f ALPHP(90a) , AbP5.‘n90 
i ) J f ftSFC 9PU) 

0 'J 1 0 I s 1 , N 

A,FFCn=ABS{F.Ff D) ■ ■ 

J FC AEFC n .l.T.B250)G3 TO 23 
4LP;iPCn=7,65E4*EXP{-4.2£5/A:EPCI)) 

Af.PRivCn=l .5f;44»FjXPC-2.3e5/AEf (D) 

3C1)=PCI )»MiJPCl)4ALPHPCIl^‘EFCi:)+ECi)*='10^CI)*AL.PHNCI)*SFCI3 

G'J TD 10 

GCI)=0.O 

coNriNJE - 

RETOKM 

£ND 

THIS SOBROUTINE DET!g,R^IMEB CDEFFIC 1 SiS f' HUTRIX FOR VARIABLES PrtlP 
SHBROliri,ME COEFF2CN,H2,FA.,S,MUP,SI, R1 , R2,3.r8,C,D,Q3 
OIMSWSI.OM MUPCN) ,SI(N) ,R1 (N) ,R2CM) , 3.C ^3, 3(M) ,C(n5 ,DCM) 

INTEGER Z . , ’ ■ 

DO 10 1=1, N 

£FCr..E3,l) GO TO 20 
IFCl.EO.N) GO TO 30 

8(1) =-:4(]ptl3*£XPC-C SI (I) +51(1-13 3/2 .0 

D(I)=-:‘^UP(I3*gXP(-CSI(I)+SI(I+l))/2'.0 

C(I)=-B(I3-D(I3 

IFU.E0.2)GO TO 50 

IFU.e2,43 GO TO 50 , 

Q(I3=H2nRlCI3+R2(I3-G(I3 3 

GO ro 10 

0(I3=H2»CRI(I3-G(I3) 

GO ro 10 

Sicr.-Us22. 263711 

GO ro 40 

SICI+1)=9.4980224+VA 

GO ro 40 

CONFIMOE 

Q(ij=9(i!-aa 

END 

THIS S0SROUTINE DETER4INES^:DEFFICIS;>ir' SATRIX FOR VARIABLSl PHIN 
SUBROUriNE C0EFF3(N,H2,VA:,Z,MUN,SI, R1, R2, G,S,C,D,Q3 
DIMSNSIOfI MUNCN) ,SI(fl) ,Rt (N) ,R2CH3 , SI M 3 , 3 ( N3 , C ( N3 , 6 C N 3 , 3 ( 4 3 
INTEGER' Z 
DO 10 I=l » N 
IFCL.E,3 .i5 go to 2(3 ' 


:n3-dcn j*EXPCtA) 



'13 


5 3 

i D 


i!) 

20 

39 

19 


13 
2 9 


ip’(r.e;3. \i) go to 30 

da)=-wia.<,ujnvxPC(SHi)+si(i-i))/2.') 

iFCa,e;5.2 )GU TO 50 
iP'Ci.SJ.'DGD TO 50 
vi(i)=-t^2^'CHl Cl) + f<2C n-:GCI ) ) 

G 3 r 0 ID 

v(I)=rt2T(Ri (D-gcD) 

GO ro 10 

Sl(r-U=22, 263711 

GO ro 40 

SUt,+ l) = 9.4980224 + VA 

G 0 r Q 4 9 

CONFIMUE 

5a) = acn-B(i) 

3C.O=5(Ni)-D(N5*eXPC"VA) 

ftErjRt^ 

Si'MD 

THIS SUBROUTINE DETSRmMSS GiURRENT OE^SITTES 
SUBROUTINE CRNTS ( N , HI , VA, MUP , 'lUN , PR?,. , P, £:, JN , JP) 

DIMENSION MUP(N),MUN(N),PHP(N) ,PHN( ,\() ,?,C!<),S,(N) JP(N: 

03 10 i:=l,N 

ifc|.,e5,i5 go to 20 

; f ^ • 


IFCt, 


N) GO TO 30 

=”;'^ypu)*PCT)tcp^p )<;c 2v?Hy 


jmci.)=‘-mum(i)*e(I)’»'(phnci:+i)-ph^^(I-ii)) /:C2,*tiU 

GO TO 10 

PHP(I-13=0.0 

PRN(I-1)=0,Q 

GO TO 10 

pHPCi+n=vA 

PHN(H-l)=VA 

GO TO 40 

CONTINUE 

RETURN 

END 

THIS SUBROUTINE DETERMINES CURRENT GAIN 
SUBROUTINE GAINCN , JP, JN . R1 , R2 , X, HFEO , 
DIMENSION JP(N3,JNCN),RlCN3rP2CM),XCN) ,RC900) 
READ' HFE,JC,JB 
JC=JPCM3+JN(N) 

DO 10 1*1 fN 

RCI)=R1(i5+R2CI3 

'^ONTTNUF 

SaDU' DOfGAFCX.R,N,JB,ER,IFAIL) 
IFCIFAIL)40,46,20 
IFCIFAIL.EQ.l) ^RITEC30,*)IFAIL 
IFaFAli,BQ.23 WRITE(30,*)IFAIL 


P ^ G£ ” 


, JN(N) 



P A 


4 D 


1 J 


ID 


2D 


iFClFf\IL,ey.3) ^P.ITEC30,*) IFAIL^ 

Si'O? 


ftSTJHN 
E ’J D 

SJBRDUTIME DETSRMINS5 .v(AX :HANGES IN POrS.NTI4LS 

SJBRD Jr^^^E DLMCN,SI,DEL,D£l.^X,D£l,MM) 

SICN),DEL(N) 

a5;i..M.^i=D.D 

DD ID 1 = 1, N 

A.DEL,«;\3SCDELCI)) 

i)E(j4K = AMAXa tDEL^X,ADEL) 

DFL'^N = 0ELMN + ADEL 

DEL^M = aE£jMN/W 

RET'JRM ' ■ 

END 

THIS SUBROUTINE DETERMINES X BY CROJT'-S FACTORISATIDfl METiOD' 

SUaROUXlNE CROUTCN,B,::,0,Q,X) 

DIMENSION 8(N) ,CCm 5 ,DCN) fOCN) ,X(N3 , IM C 9 3 D ) , PM t 900) 
HM(U=-0(U/C£i) 

PM(l)=3C13/CCl) 

DO 10 !=2rN 
HMr“' ' ■ 

PM( ^ 

CONTINUE 
X(N)=Py|(N) 

N1=N-1 

00 20 1=1, N1 
J=N-I 

KCJ)=PM( J3+HMCJ)*XC J+i) 

CONTINUE 
RETURN 
END 


i 10 1=2, N _ 

l(I)=-0(i)/CC(I)+BCI)=t=HM'CI-U) _ , 

tc I) = C 0(1 3 -8( I l^PMC I-;l ) ) / C CiC I ) + B( I ) << IMC I-l 3 3 



I- t f- i * t f .f .f. ^ 4 '4? ^ f f:- # ■ f- 


■ . ■ 3 M. /!•• :W’^S 1 

:4; I 4£ 4 5f.' 4' |: ^ 


U 

23 

33 

43 

53 

:i 

73 


133 


PRDCaA'»> M3. 2 

THIS PROGRAM DETERMIMES SMiriER SGRPAGE: GOMGENTRATI L);'i S 
REaJxRSD FOR EQUAL SMITTSiR GHARGE F3R OIFFEREMT EMITTER 
Ve,R}h(iLF, OGGLARATTDM . 

REAL n'.LprNl,HUXE,NS.,M3,NG:f BETA, ALi?RA,)CC10U) ,DaP(100D , 

lf5f(iD3).Q£,EK 

iUTSGF.R Mf],fFAIL 

GDNSIAMrS DECLARATION 

VT=3,326 

L(D=3'4,E-4 

MI = !.5S10 

;«,3F POINTS COt^SIOERE-D FOR IMTEGRATIO.^ 

MNsRP 

OUTPUT FILE 

3PE'i(UMir = 30,DEVIGE='DSK',FrLS='QE.'l ') 

CONSIDERATION OF DIFFSREJIT SMiTrER DEP, T-IS 
FOR IMUM EMITTER 
K=2 FOR 2MUM EMITTER 
K=3 FOR 4MUM EMITTER 
iC=4 FOR 4.75MUM EMITTER 
K=5 FOR 8MUM EMITTER 
03 233 K=2,2 
IFC<„E3,UG0 TO 10 
23GO TO 20 
3)G0 TO 30 
43GO TO 40 
,.5)00 TO 50 
!-4)/L0?G0 to 70 
KE=C2.E,-4)/LO;GO TO 70 
XEs:C4.E-4)/LO|GO TO 70 
XE=C4.75E-43/LD?G0 TO 70 
XE=C8.S-4}/LD?GO TO 70 
DOPING PROFILE PARAMETERS 
ME* (7. El 9) /Ml 
MB=C2.E1?5/N1 
MC=C2.S145/NI 
8ETAsl '2 

ALP4A = AL0G(NE/(N8tEXP(-8E.TA)-MC)) 

STEP DISTRIBUTION 

hi = xe/:mm 

X(1)=0.0 
DO 100 1=1, NN 
J = Ifl 
x(j) = r*'Hi 
CONflMUE ■ ■ 


ifck.e: 
iFc K. e; 
if(k.e: 
ifcKje: 
XE=Ci 


DEPTHS 

EMTC1003 





i?.) 


1 %’) 


•ij 

20 

u 


>4 s ,'J ’'a +■ 1 

aSXr^aM^,M4T^^^4 OF DOPING PROFILE 

L'tDP(.r) = 'ji=^fexP(-ALPHi!k<'(CKCl)/X£3*^2) )-:M3*£XPC-BErft^C (KCI j/iCS)^»2 J 

1 ) t'3C 

CDNflMlJE 

.•'IIFFFCN, VT,NI,DOP,S-HJT) 

DO 320 I. = lfN 
x(i)=DOP(i)/eNi(i) 

COi^l’i'JLJF 
iFftri,=i 


:/\L0 001GS4F(X,y .N.QE,£R, IFAIL) 
IF(r,FML) 440 .<l 46,330 
^iRirE(30,*)IFAIL 


XE=XE*LDJ ■•'iE = NE*NI;«RirE(30,3')' KE, NS,, 3'Si, ER' 

CONriM'JE 

SfQ? 

END ■ ' 

rniS S08R0UTINE DErE.R4INS5 3GM OORREGinON TO NI 
SUBROUriNE NIEFFCN,/T,NI,OOP,ENi:) , 

DIMENSION OOP(N),ENnN) ,FCiD6) fADOPClOO) 

REAL NI 
00 10 Isi,N 

XL=C1 .E.175/NI 
ADOP(I 3 sA 8 S(DOP(I )3 
IFCAOOPd) ,Ge.XL3 GO TO 20 
ENIC,I3 = 1. 

CONTINUE 

RETURN 

END ■ 



